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\ We study a spin chain with an anisotropic XXZ coupling in a magnetic field. Such 
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a chain models several proposed types of quantum computers. The chain contains a 
defect with a different on-site energy. The interaction between excitations leads to 
the onset of two-excitation states localized next to the defect. In a resonant situation 
scattering of excitations on each other might cause decay of an excitation localized on 
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J> ■ the defect. We find that destructive many-particle quantum interference eliminates 

\D , 

. this decay. Numerical results confirm the analytical predictions. 
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■ An important motivation for making a quantum computer as envisaged by Feynman is 

O ■ to simulate real quantum systems. Disordered systems with interacting excitations are of 
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particular interest in this context. An important class of such systems are disordered spin- 



^ ' 1/2 systems. The spin-spin interaction leads to a complicated excitation spectrum, that has 



been described only in special cases where the system is integrable 

A spin-1/2 system can be modelled by a set of interacting two-level subsystems, qubits. 
Controlled qubits, with interqubit coupling that is not turned on and ofL are basic elements 
of several proposed implementations of a quantum computer (QC) . Qubit 

energies can be often individually varied. This allows one to explore resonant situations 
in a multiply-excited system. In contrast to naturally occurring many-particle resonances 
in condensed-matter systems, in a QC the parameters can be tuned into resonance in a 
controlled way. 
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In this Letter we show that a simple spin chain with a defect displays new types of 
localized two-excitation states and destructive two-excitation interference, i.e., an effect of 
two-excitation antiresonance. This is a result of the interplay between disorder and inter- 
action. In turn, localization of excitations on individual qubits is important for quantum 
computing, because the qubits can be then prepared in desired states, which are not de- 
stroyed by other excitations. 

We consider a spin chain in a magnetic field. The defect corresponds to "diagonal dis- 
order": the spin energy on the defect site differs from that on the host site (for example, 
the magnetic field on the defect site differs from the field on the host sites). The exchange 
coupling is anisotropic and is described by the XXZ model, which applies to several proposed 
quantum computers. The Hamiltonian of the chain with a defect on site no is of the form 

if = ly e^a^ + ly y Jii<K+^, (i) 

£ ~\~ 9Sn,noi Jxx Jyy Jzz t/A. 

Here, are the Pauli matrices and e*^"^ is a spin-flip energy on site n. All e^'^'^ are the same 
except for the site n = no- This is not the condition in which a QC should normally operate, 
because excitations are mostly delocalized and are scattered by each other. But even here, 
as we show, a multiply-excited system has strongly localized states. 

In Eq. (^, the defect is characterized by the excess energy g . The parameter A deter- 
mines the coupling anisotropy. We assume that |A| ^ 1, as in the case of a QC based on 
electrons on helium, for example j^. To simplify notations, we further assume that the spin 
coupling is antiferromagnetic: J, A > 0. 

For large on-site energies e^^\ the ground state of the system corresponds to all spins 
pointing downward, (cr^) = — 1. The number of excitations is the number of spins pointing 
upward. It is conserved, because the total spin projection J2 <^n commutes with H. In what 
follows we count energy off from the ground-state energy. 

Single-spin excitations for the Hamiltonian are well understood j^. In an infinite 
chain they are either of the spin-wave type ("magnons") or localized on the defect. The 
magnon energies form a band of width 2 J centered at £i = e — J A. The localized excitation 
has energy ei + {g'^ + J'^Y^'^ sgng. We will be interested in the case g ^ J, when the excitation 
is strongly localized on the defect. 

Two-spin excitations are well understood for an ideal chain For strong anisotropy, 



3 



A ^ 1, they are either uncoupled magnons or bound pairs (BPs) of flipped neighboring 
spins that propagate together. The two-magnon energy band is centered at 2ei and has a 
width of 4J. When excitations are on neighboring sites, their energy is changed by J A ^ J, 
because of the term ( JA/4) X^n '^n'^n+i ^ O- Respectively, the BP band is centered at 
-Egp ~ 2ei + JA and is well separated from the two-magnon band. Motion of a BP occurs 
via virtual dissociation, where one excitation hops to a neighboring nonexcited site and then 
the second excitation moves next to it. The bandwidth resulting from this process can be 
calculated by perturbation theory in 1/A. It is given by J/A and is parametrically smaller 
than the magnon bandwidth. 

In the presence of a defect, there emerges an additional two-excitation state where one 
excitation is localized on the defect and the other is in an extended (magnon-type) state. 
We call it a localized-delocalized pair (LDP). The LDP band is centered at -Eldp ~ + g 
and has a width of the one-magnon band 2 J (see below). 

If the LDP energy band overlaps with the BP band, one would expect that the two types 
of excitations mix together. In other words, the propagating magnon of an LDP can scatter 
off the localized excitation, and they together would move away from the defect as a BP. 
This would lead to magnon- induced decay of a localized excitation. We show below that the 
decay does not happen because of destructive quantum interference. 

The wave function of a chain with two excitations is a linear superposition ^2 = 
Z]n<m '^('^j ^)) where (f){n,m) is the state where the spins on sites n and m are 
pointing upward and all other spins are pointing downward. From Eq. (^, the Schrodinger 
equation for the coefficients a{n, m) has the form 

+ JA6m,n+i) a{n,m) 
+ ( J/2) J2k=±^^^^ + k,m)+ a{n, m + k)] = E2a{n, m). (2) 

Here, E2 is the energy of a two-excitation state counted off from the ground-state energy, 
and a{n,n) = (two excitations may not be placed on one site). 

The defect may lead to a localized state of a bound pair with energy close to the energy 
of a propagating pair E^p. This state is centered on sites [uq + 1,^0 + 2), or equivalently 
(no — 1, no — 2). To show this we seek the major term in the solution of Eq. Q in the form 

aW(n,m) = A6n,no+iSm,no+2 + Cie^^^("~"°-2)5„,„+ie(n - % - 2) 

+ (c'2e'^2(™-"°-2) + c^e-^^'^^""""-')) 5n,noQ{m - uq - 2). (3) 
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Here, A is the amplitude of the pair on sites (no + 1, no + 2), Ci is the amphtude of a BP 
propagating away from the defect with wave number 9i, and C2 and are the amphtudes 
of LDP waves in which one excitation is on the defect whereas the other is propagating away 
from and toward the defect, respectively, with wave number 62; 0(m) is the step function, 
it is equal to 1 for m > and for m < 0. For concreteness, we consider excitations "to 
the right" from the defect, m > uq. They do not mix with the excitations "to the left", to 
first order in A^^. 

The BP localized next to the defect would be described by the wave function Q with 
6*2 = and with Im 61^2 > 0. On the other hand, resonant mixing of LDPs and BPs would 
be described as scattering of an LDP wave incident on the defect (Cg = 1,62 > 0) into a 
reflected LDP with amplitude C2 and a BP with amplitude Ci and 61 > 0. 

Eqs. (121), © with n = m — l>no + 2 and n = no, m > uq + 3 give, respectively, the 
dispersion laws of BPs and LDPs far from the defect. From Eq. (0), a BP on sites (n, n + 1) 
is directly coupled to nonresonant dissociated pairs on sites (n — l,n + 1) and (n, n + 2). 
Similarly, an LDP on sites (no,m) is coupled to resonant pairs on sites (no,m ± 1) and to 
nonresonant pairs on sites (no ± l,m). Finding the amplitudes of the nonresonant pairs to 
first order in A"^ from Eq. ^ and substituting them back into the equations for Ci and 
C2, we obtain, respectively, the BP and LDP dispersion laws: 

EbAO) = ^BP + (^/2A) cos 9, = 2^1 + + ( J/2A), (4) 

ELDp{e) = 4dp + ^cos^, 4dp = 2ei + g+ (J/2 A) 

(the renormalization of the energy spectrum is calculated for \E2 — E^p\ ^ J). 

Wave mixing near the defect can be described in a similar way by writing Eqs. (0), 
for the sites (no + 2, no + 3), (no + 1, no + 2), and (no, no + 2) and eliminating nonresonant 
pair states by perturbation theory. For energies E2 close to E^p this gives 

Cie-'^' =A + C2e'^' + C^e"'^^ , A = Cze"'^^ + C^e'^^ + A"^ (C2 + C^) , (5) 

4(JA -g- Jcos02) - JA"^] A + 2J{C2 + C2) 

+JA-i(C2e^^2 + C'^e~'^^) + JA-^Ci = (6) 

(the corrections ~ A^^ to the first equation (jSJ come from higher order terms of the per- 
turbation theory; they are not important for the present calculation). 
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Eqs. (0), © contain two dimensionless parameters, A and a = {J A — g)/J, which are 
related to the three energy scales in the system: the distance between the BP and LDP 
bands J A — g and the LDP and BP bandwidths J and J/ A. 

To find the state localized on sites (no+1, no+2), we set C2 = 0. Then, for A ^ 1 the solu- 
tion of Eqs. (jni), & can be approximated by exp(—-i6'2) ~ 2J^^(JA — (7) + (2A)~-'^[exp(26'i) — 1] 
(the term oc A~^ is important only for |a| ^ 1 and should be modified for |a| ~ 1). The 
energy of the localized state i^ioc and the reciprocal localization lengths Im 612 are obtained 
by substituting this expression into the equation Eioc = -E'bp(6'i) = -Eldp(6'2)- 

For |a| ^ 1 the BP and LDP bands are far from each other, and 

The localized state amplitude is A. The state is mostly hybridized with the BP band, with 
Ci/A a/{A — a) and Im 9i = ln[(A — a)/a]. The admixture of the LDP band is small, 
~ l/2a; ^ 1. In a way, the localized state is a "surface" state of the BP band, with 
the role of the surface being played by the defect. It emerges for a < A/2. 

As the BP and LDP bands approach each other and a decreases, the energy Ei^c moves 
away from the BP band. The localized state becomes stronger hybridized with the LDP 
band than with the BP band, with C2/A ^ l/2a and Ci/A ~ a/A < 1 for |a| < A. The 
energy separation from the LDP band is 

^loe-4DP~^[«+(4«)"'] (8) 

It is of the order of the LDP bandwidth J. The localized state exists for |a| > 1/2. The 
reciprocal localization length Im 62 = ln(2|Q;|) goes to zero as \a\ —>■ 1/2. As the BP band 
goes through the LDP band with varying J A — g, the localized state first merges with the 
LDP band and disappears for a = 1/2, and then is split off from the opposite side of the 
LDP band for a = —1/2. In more details the dependence of i^ioc on parameters is discussed 



elsewhere 







along with the onset of other localized two-excitation states. 
We now consider the possibility of resonant mixing of LDP and BP states where the 
bands overlap. As mentioned above, we can do it by studying scattering of an LDP magnon 
off the excitation on the defect, i.e., by setting = 1 in Eqs. ©, (©• 

The LDP band is much broader than the BP band. From Eqs. (@]), LDPs with energies 
within the BP band have wave numbers 62 such that cos 6*2 ~ {J A — g)/J. For such LDPs, 
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C2 ~ — C2 = —1 from Eq. i.e., they are fully reflected from the defect. The amplitude 
of the BP wave is Ci = 0, to zeroth order in A~^. This happens even though the BP on 
sites (no + 1, no + 2) has a large amplitude A ^ 2i sin6'2. 

The vanishing of the BP wave is a result of destructive interference, or antiresonance, as 
seen from the first of equations The BP wave is coupled to both the pair (no + 1, no + 2) 
and the pair on sites (no, no + 3), which has an amplitude C2exp{i62) + C2 exp(— 'i6'2). For 
resonant values of 62 the complex amplitudes of these two pairs are equal and opposite in 
sign, they cancel each other. One can show in a similar way that if a BP wave is incident 
on a defect, it will be fully reflected and no LDP waves will be excited. Absence of mixing 
of BP and LDP states at resonance is illustrated in Fig. which is obtained by direct 
diagonalization of Eq. (j21). 




FIG. 1: Amplitudes of two-excitation states on sites (n, m) in a periodic chain. Left and right 
panels show, respectively, bound pairs and localized-delocalized pairs with nearly same energy. 
The chain length is 20, the defect is on site no = 5, and g/J = A = 10. The BP is located on 
sites m = n + 1 and is well approximated by a{n,n + 1) oc sin[107r(n — 6)/17]. Because of the 
antiresonance, the BP amplitude is equal to zero on sites no, no =b 1, no ± 2. The LDP is located on 
sites {nQ,m), with m — uq > l(mod20), and is well described by a{nQ,m) oc sin[(m — 1)tt/2]. The 
LDP is strongly hybridized with the pairs on sites (no + /c,no + 2k) (k = ±1), which have same 
amplitude, as expected. 

In this paper the effect of antiresonance has been found for many-particle excitations. 
It leads to localization of an excitation on a defect in an XXZ spin or qubit chain even 
where the excitation would be expected to decay via multiparticle scattering. We also found 
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new localized two-excitation states. Their energy and localization length depend on the 
interrelation between the interaction anisotropy A^^ and the relative excess energy of an 
excitation on the defect site g/J. The results indicate that quantum computers can be 
advantageous for studying new many-body effects in systems with disorder. 

From the standpoint of quantum computing, the many-particle localization that we dis- 
cuss is important, because it is much easier to manipulate excitations when they are localized 
on individual qubits. Of interest for quantum computing is also the occurrence of different 
types of localized two-excitation states on neighboring qubits. It allows creating on these 
qubits entangled Bell-type states, which are to some extent similar to recently studied en- 
tangled excitonic states in a quantum dot [12], but which appear in the present case in an 
extended system. We expect that coherent multi-qubit states can be created in a chain with 
a defect as well. 
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